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ABSTRACT 

We cross match the NVSS and FIRST surveys with three large photometric cat- 
alogues of luminous red galaxies (LRGs) to define radio-loud samples. These have 
median redshifts 0.35, 0.55 and 0.68 and, by matching rest-frame optical and radio 
properties, we construct uniform samples across the three surveys. This paper is con- 
cerned with the clustering properties of these samples derived from the angular cor- 
relation function. The primary aim is to characterise any evolution in the clustering 
amplitude of radio galaxies bellow z ~ 0.68. 

We find no evidence for evolution in the large-scale (~ 1 — 50 hr 1 Mpc) clustering 
amplitude. Our radio galaxy autocorrelations are consistent with previous findings 
indicating little-to-no evolution in the redshift range 0.68 to (~ 6Gyr of time). We 
also cross correlate radio galaxies with the parent LRG samples to increase the preci- 
sion of our results and again find no evidence for evolution in comoving coordinates. 
Our results are inconsistent with a long-lived model for the clustering evolution that 
assumes radio sources randomly sample the LRG population. A model where the halo 
mass is constant with redshift is consistent with the data. This is similar to QSOs 
that have clustering amplitudes consistent with a single halo mass at all redshifts. 
Given that the brightest radio sources show stronger evolution in space density com- 
pared to fainter radio sources we restrict our samples to include only objects with 
L > 10 26 W/Hz and repeat the analysis. Again we find no evidence for evolution in 
the comoving correlation amplitude. These radio sources appear to inhabit the same 
mass halos as fainter radio galaxies (~9x 10 13 /i _1 M©). These halos are ~ twice as 
massive as those of the general LRG population (~4x 10 13 /i -1 M©) and ~ 30 times 
as massive as optical AGN/QSOs (~3x 10 12 /i _1 M©). 
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1 INTRODUCTION 

Radio galaxies have long been used to study large scale 
structure as they are highly biased indicators of the mass 
distribution and are naturally found in the high redshift 
(z ~ 1) Universe. Understanding the observed bias of ra- 
dio galaxies requires relating their clustering to that of their 
host populations. Furthermore, any differential evolution in 
their clustering with respect to their hosts gives insight into 
the nature of radio galaxies and the mechanisms that trigger 
them. 

The typical host galaxy of a radio source depends 
strongly on its luminosity. At lower radio luminosities (L < 
10 23 W/Hz) the population is almost exclusively made up 
of starforming galaxies with related diffuse radio emis- 
sion, above L ~ 10 23 W/Hz almost all radio sources are 
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associated with active galactic nuclei (AGN). The hosts 
of th e brightest radio galaxies have strong emission lines 
(e.g. iHine fc Longaiil 1 19791 : iMcCarthvl 1 19931 ) while lower- 
luminosity radio AGN are hosted by massive ellipticals 
with little-to-no optical line emission. This transition im- 
plies a cutoff luminosity between object s with and with- 
out emission lines that Ijohnston et al.l (J2008) found to 
be around ~ 10 26 W/Hz. This luminosity is close to the 
traditional cutoff b etween morphologically classified FRI 
and FRII objects (|Fanaroff fc Rilevl 1 19741 ). Furthermore, 
while the bright end of the radio luminosity function 
evolves strongly towards higher space densities at high z 
|Longairlll96llDunlop fc Peacocklll990l ). fainter radio AGN 
L 5, 10 26 W/Hz) show much less eviden ce for evolution 
Clewlev fc Jarvisj|2004l : ISadler et al.ll2007l ). This evolution 
in space density is roughly matched by the typical hosts of 
radio AGN. Massive elliptical galaxies (luminous red galax- 
ies; LRGs), that host lower luminosity radio AGN, show lit- 
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tie evolution in number density for z < 1 IjBell et al.l 12004 ; 
IWake et all 120061 ; [Brown et alj|2007h . Higher luminosity ra- 
dio sources are more associated with optical/Xray AGN 

McCarthy 1993) th at show strong evolution out toz~ 2.5 

Schmidtlll968l : ICroom et al.ll2009l ). 



Radio sources are strongly biased tracers of the un- 
derly in g mass distribution o f the Universe II Yates et al] 



19891: iPeacock fc Nicholson! Il99ll; 
Magliocchetti ct al. 2004; Brand ct al 



Blake fc Wall! 
20051 ) 



2002; 



Since the ad- 
vent of deep all sky surveys (e.g. FIRST, NVSS, SUMSS, 
WENSS) there have been major leaps in the ability to de- 
scribe the environments of radio galaxies. These milli-Jansky 
surveys sample the luminosity range ~ 10 23 to 10 26 W/Hz 
at moderate redshifts, where radio sources are found almost 
exclusively in LRGs. LRGs are known to be a strongly bi- 
ased tracer of the matter density, and their bias increases 
with the mass/luminosity of the LRG. 

Radio sources tend to be found in the most massive el- 
liptical galaxi es and more massive galaxies are more likely 
to host them JBest et al.ll2005l ; iMauch fc Sadlerll2007tl . Re- 
cent studies have compared the clustering of radio galaxies 
to that of radio-quiet galaxies that have been matched in 
their optical properties to the radio sample. Most find that 
radio galaxies are significantly more clust ered than optically 
identical samples of quiescent galaxies (|Wake et al.l 120081 : 
Mandelbaum et al.ll2009l ; iDonoso et all 120091 ') . although see 
Hickox et al.l (]20091 ) for a counter example we will discuss 



later. The indication is that the environment of a galaxy 
contributes to the probability of it hosting a radio source. 

In this paper we are primarily concerned with how the 
clustering of radio sources evolves for z < 0.7. In this redshift 
interval t he clustering ampl i tude of LRGs is approx imately 
constant l|Wake et al.l 120081 ; ISawangwit et alj|2009h . How- 
ever, the clustering amplitude of quasars evolves such that, 
when models of gravitational collapse are assumed, the im- 

J 'lied dark halo ma s s for quasars is a pproximately constant 
Croom et al"1l2005l : IRoss et alJl2009l ). We will evaluate the 
clustering amplitude of radio LRGs in a series of samples 
from z ~ 0.68 to 0.35 to investigate any evolution in their 
clustering. This will be done both for a sample of medium 
luminosity (L ~ 10 24 ' 7 W/Hz) radio galaxies, and a high 
luminosity subsample (L > 10 26 W/Hz). 

In section [2] we introduce the data that will be used 
in this work, section [3] describes the algorithm we develop 
for identifying radio LRGs, section [4] describes the corre- 
lation analysis to be used, section [5] compares the cluster- 
ing of radio-loud and quiet LRGs, section [6] then looks at 
the evolution of radio LRGs and in section [7] we investigate 
the evolution of the very brightest radio LRGs. Through- 
out this paper we assume a flat (fi m ,f2A) = (0.3,0.7), 
Ho = 70kms -1 Mpc~ 1 cosmology. All radio flux densities 
and luminosities are at 1.4 GHz unless otherwise stated and 
when estimating radio luminosities throughout this paper 
fc-corrections have been performed assuming a continuum 
shape of S v oc v~ ' 7 . 



2 DATA 

The LRG samples used in this paper are originally de- 
fined from three spectroscopic surveys: SPSS, 2SLAQ and 
A An (|Eisenstein et al.ll200ll : ICannon et~aLll2006l ; lRoss et all 



2008). The LRG selection was refined by ISawangwit et al.l 
1 2009 ) to create photometric samples from the SP SS DR5 
( Yorket al J 12000 : lAdelman-McCarthv et"aD 120071 ) and to 
cut down on stellar contamination. We use LRG samples 
defined identically but drawn from the more recent PR7 
IJAbazaiian et al-lfeoosl ). 

To identify radio-loud LRGs in our sample we com- 
pare our LRG catalogues with radio source catalogues from 
the N VSS (|Condonet al J 1 19981 ) and FIRST (JBecker et all 
120001 ) surveys. Both surveys are carried out at 1.4 GHz and 
the high angular resolution (~5 arcsec) and faint flux limit 
(la ~0.15mJy) of the FIRST survey combined with the 
large-scale sensitivity of the NVSS make them complemen- 
tary tools for identifying radio LRGs. The FIRST survey 
covers much of the same sky area as the SDSS north Galac- 
tic cap. In the overlap region there are 110104, 652401 and 
799519 objects in t he SDSS, 2SLAQ andAAfi LRG sam- 
ples respectively. See lSawangwit et al.l |2009) for the redshift 
distributions, expected contamination and autocorrelation 
clustering properties of the LRG samples. 



3 RADIO CROSS MATCHING 

Complex radio source morphology due to extended struc- 
tures compounded by size-dependent response effects due 
to the use of interferometers makes radio source matching 
more involved than simply matching sky positions. The most 
accurate and precise method for cross-matching has always 
been manual inspection (e.g. ISadler et al.l 120071 ). However, 
the size of our LRG samples forces us to define an auto- 
mated cross-matching procedure that does not require visual 
inspection of all of the potential radio matches in the sam- 
ple. Automated radio-matching routines have been deve l- 
oped before (e.g. iBest et"afl 120051 : [kimball fc Ivezidl2008l ). 
and even samples defined by manual inspection use some au- 
tomated procedures to define a sample of potential matches 
that are then visually inspected. 



3.1 Our cross-matching procedure 

ISadler et al.l l|2007l ) define two samples of radio LRGs from 
the 2SLAQ spectroscopic sample. One is based on manual 
inspection and the other based on automatic cross matching 
to the FIRST survey. We modified their selection criteria to 
define a radio matching procedure and use their manually 
inspected sample to test our results. We match the LRG and 
FIRST catalogues within a 30 arcsec radius, and the LRG 
and NVSS catalogues within a 180 arcsec radius, retaining 
all radio matches within the given radius around an LRG. 
In this initial cross matching we do not app ly any radio 
flux li mits, these are applied below. Following I Sadler et al.l 
(2007) we use the FIRST matches as our primary tool for 
identifying radio galaxies. We accept matches that meet one 
of three criteria: 

i) A FIRST match is within 3 arcsec of the LRG position. 
ii) A FIRST match is within 10 arcsec of the LRG position 
and the major axis of the source is orientated within 25° of 
the angle to the LRG position. 

ill) Given more than one radio match, the flux-weighted 
mean of the FIRST positions (or any subset thereof) is 
within 6 arcsec of the LRG position. 
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For any sources that had NVSS matches within 
180 arcsec but not FIRST mat c hes w e applied a similar set 
of criteria based on iBest et al.l (J2005|) to the NVSS sources 
to define potential radio galaxies. These are: 

i) An NVSS match is within 10 arcsec of the LRG position. 
it) A galaxy matched to two radio components that is 
< 10 arcsec from the flux weighted radio centroid, and the 
radio-galaxy-radio angle is > 150° 

Hi) Given more than two radio matches, the flux-weighted 
mean of the NVSS positions is within 10 arcsec of the LRG 
position. 

We apply this procedure to the 2SLAQ spectroscopic 
catalogue and find 428 objects satisfy our selection criteria. 
Of these, 57 objects do not appear in the Sadler catalogue, 
while 20 objects from the Sadler catalogue are not found. 
This indicates that we are ~ 5 % incomplete with respect to 
the Sadler catalogue, the majority of sources missed by our 
selection criteria have complex structures in their FIRST 
images. Of the 57 extra matches 29 can be explained as ob- 
jects that had NVSS matches but no corresponding FIRST 
sources. These objects were rejected by Sadler et al. to im- 
prove the reliabil ity of their c atalog ue. We take an approach 
more similar to IBest et al.l 1)20051 ) including these sources 
given the criteria above. This indicates that our routine, ap- 
plied to the 2SLAQ spectroscopic LRG catalogue, is ~ 93 % 
reliable. 

We apply this radio cross-matching routine to the 
SDSS, 2SLAQ and AAf2 photometric LRG samples and 
obtain 9689, 15236 and 8364 radio LRGs from the three 
samples. This corresponds to 8.8, 2.3 and 1.0 % of the par- 
ent samples. The considerably higher detection rate for the 
SDSS sample is due to the lower redshift and higher absolute 
(optical) magnitude of the SDSS LRGs 

We perform a simple check of the reliability of our final 
radio samples by offsetting the original LRG samples, then 
rerunning the cross-matching code on the offset catalogues. 
We offset the samples by ±0.1° in declination (twice the 
180 arcsec maximum cross-matching radius used). The cross- 
matching code found 171/138, 516/453, 457/499 matches 
in the offset SDSS, 2SLAQ and AAfi samples respectively 
for the positive and negative offsets. This corresponds to 
~1.5-6 % contamination in our final cat alogues similar to th e 
~ 7% value found when comparing to lSadler et al.l (|2007l ). 

Many of the radio identifications are based solely on 
matching to FIRST source catalogues thanks to the angu- 
lar resolution of the FIRST survey. However, the FIRST 
survey's smaller beam resolves out large-scale structure and 
does not give accurate radio powers for these sources. We 
cross match our final radio LRG catalogue with the NVSS 
within a 180 arcsec radius. We then identify which NVSS 
matches are associated with the galaxy using the criteria 
above. We sum these NVSS components to give the radio 
flux densities used throughout this pa per. The NVSS is es - 
sentially complete above ~ 2.8 mjy II Condon et al.l Il998l ) 
and in our final sample we apply a radio flux density cut 
of 3 mjy. 

3.2 Redshift distributions 

The radio-loud subsets of the three LRG samples may not 
have the same redshift distribution as their parent sample. In 
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Figure 1. The solid lines shows the redshift distributions of the 
SDSS, 2SLAQ and AAf! spectroscopic samples (left to right). 
The dashed lines show the redshift distribution of those objects 
identified as radio sources. These are, in general, identical to the 
parent distribution. In the highest redshift AAQ sample there is a 
difference, but there are only 17 objects defining this distribution. 



Fig. [T] we plot the redshift distribution of the spectroscopic 
LRG samples, and those objects that we identify as radio 
galaxies from our matching routine. 

For the SDSS and 2SLAQ samples it is clear that the 
radio-detected subsamples have an identical redshift dis- 
tribution to their parent samples. The AAf2 radio sample 
shows indications of being at a higher redshift. However, 
there are only 17 spectroscopic AAfi radio LRGs and hence 
the distribution is highly uncertain. For the rest of this paper 
we will assume that the radio LRGs have the same redshift 
distribution as their parent sample. None the less, if the in- 
dication of Fig. [l] is accurate (that the radio LRGs have a 
higher redshift than the parent sample) this will effect our 
results. In particular a cross correlation between the radio 
loud sample and the parent sample will be reduced if the 
redshift distributions do not overlap so well. 



4 CORRELATION ANALYSIS 

To study the clustering of the samples defined above we cal- 
culate the 2-point angular correlation function to measure 
the LRG x LRG autocorrelations and the LRG x radio cross 
correlation. The 2-point angular correlation function is de- 
fined as the relative increase in pair counts at a given scale 
over that of a random sample. That is 



_ N(R) DD(6) 
W{0) ~ N(D) DR{6) ~ 1 



(1) 



where DD gives the number of sample-sample pairs and DR 
gives the sample-random pairs at separation 0. The expres- 
sion is normalised by the number density of randoms N(R) 
divided by the density of the sample N(D). In an autocor- 
relation the data-data pairs come from the same sample, in 
a cross correlation t hey come from differing samples. Othe r 
estimators for w(0) l|Hamiltonlll993l ; lLandv fc Szalavlll993l ) 
give equivalent results for this analysis. 

We constructed our random catalogue from the DR7 
masks published on the SDSS website. The full DR7 area 
was then cut down to the area covered by FIRST and popu- 
lated with random sources with a density twenty times that 
of the radio LRGs. 
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4.1 Error estimation 

Standard Poisson errors are notoriously inaccurate for cor- 
relation f unctions due to the high level of covariance in the 
data (e.g. ISawangwit et al.ll2009T ). More often bootstrap or 
field-to-field errors are employed as an independent measure 
of the precision of the results. We calculate field-to-field er- 
rors by subdividing our sample into nine declination strips 
that contain even numbers of objects in our random cat- 
alogues. We calculate the correlation function within each 
field and the error is estimated by the rms of the correlation 
functions in each field divided by three. 



5 CORRELATION FUNCTIONS 

In Fig. [2] we show the LRG-LRG and LRG-radio (hereafter 
LxL and LxR) angular correlation functions for the three 
samples we are studying. In the top panels we show the raw 
correlations. In the bottom panels we show the ratio of our 
measurements to the best doub le power- law fit to the L xL 
correlation function taken from ISawangwit et al.l i|2009l ). 

In all three samples we study it is apparent that the 
LxR correlation is stronger than the LxL autocorrelation, 
indicating that radio galaxies are found preferentially in 
richer environments than their parent LRG population. The 
LxR cross correlation (like LxL) exhibits a clear inflection 
at intermediate scales (~a few arcmin) which is generally in- 
terpreted as the crossover between the intra and inter-halo 
scales. Similar to previous authors we find that the excess 
clustering signal f rom the LxR correl ation is strongest in 
the 1 halo regime. iDonoso et al.l (|2009T ) interpret this as in- 
dicating radio galaxies hold a particularly central position 
within their dark matter halo. At larger scales the difference 
between the LxL and LxR correlation functions is reduced, 
and becomes negligible above ~ 20 arcmin in the AAf2 sam- 
ple. As we discuss later this may be due to the radio sources 
being brighter in this higher-redshift sample. 

To relate the measured angular correlation function, 
w(0), to the 3D spatial correlation f unction, £(r), w e follow 
the standard approach outlined in iPhillipps et al.l II 19781 1 . 
Given the clear inflection in the correlation functions shown 
in Fig. [2] we assume the spatial correlation function can be 
parametrised as a double power la w and integrate equa- 
tion 13 from IPhillipps et al.l (|1978T 1 to give w(9). In our 
analysis we fix the break scale in the correlation function to 
1 /i _1 Mpc. We fit for the slopes of the two power laws and 
the amplitude of the large-scale power law. Table[T]gives the 
parameters of fits to the observed data. 

The large scale slope of £(r) in the LxR data is essen- 
tially identical in each sample, as it is when compared to the 
LxL correlations. The small scale power law is considerably 
steeper in the radio correlations indicating that the local 
environment of radio LRGs is more overdense compared to 
the environment external to its own dark matter halo. 



5.1 Radio loud vs. quiet LRG clustering 

Radio galaxies are known to be fo und preferentially in th e 
most luminous host galaxies (e.g. iMauch fe Sadler! |2007| ). 
The most luminous galaxies tend to also be found in the 
highest density regions and following previous authors (e.g. 



Table 1. Parameters of power-law fits to the LxL and LxR 
angul ar correlation fun ctions. The first columns are taken from 
iSawangwit et al.ll2009l) . The fits are shown in Fig. [2] 



LxL 
Sample ro(/i Mpc) 7 



LxR 

r (h~ 1 Mpc) 7 



SDSS 



7.35±0.08 
9.15±0.16 



2.19±0.03 
1.85±0.04 



7.46 
9.48±0.21 



2.31±0.07 
1.79±0.08 



2SLAQ 



6.32±0.03 
7.78±0.05 



2.16±0.01 
1.85±0.02 



6.27 
9.13±0.17 



2.32±0.04 
1.82±0.06 



AAfi 



5.96±0.03 
7.84±0.04 



2.14±0.01 
1.81±0.02 



5.06 
8.71±0.18 



2.49±0.05 
1.83±0.07 



IWake et al.ll2008l ) we test whether this covariance is respon- 
sible for the increased LxR correlation strength. We con- 
struct random samples of LRGs from our parent catalogues 
that have the same (r-band) magnitude and r — i colour dis- 
tributions as the radio LRG catalogues. In Fig. [2] we show 
the cross correlation between this matched catalogue and 
the full LRG sample. We find that they follow the LxL lines 
almost exactly indicating that the host mass/luminosity is 
not the single factor that defines the probability that a LRG 
is a radio source, there must be an environmental effect as 
well. 

Our results ar e in agreement with most previous authors 
llWake et all |200S|; [Mandelbaum et all 120091 ; IDonoso et all 
l2009h . iHickox et al.l 1J2009I '). on the other hand, finds that 
the radio galaxies in his sample cluster with the same ampli- 
tude as galaxies matched in terms of their optical properties. 
They note that the result is marginal but also point out that 
their radio galaxies have considerably higher stellar mass to 
radio power ratios than the other samples, and suggest this 
as a potential explanation for the discrepancy. Certainly if 
the clustering strength depends on the radio luminosity that 
could potentially explain their result. 



5.2 Clustering strength vs. radio luminosity 

iPrestage fc Peacockl II 19881 1 found that less luminous FR1 
objects tended to be in overdense regions while more lumi- 
nous FR2s had a clustering environment similar to the over- 
all populat ion of radio-q uiet galaxies. However, many au- 
thors (e.g. iPeacock fc Nicholson! 1 1991; Maglioc chetti et al.l 
2004) have found no relationship betw een rad i o pow er and 
large scale clustering strength. IDonoso et al.l (|2009T ) found 
variations in the correlation amplitude with radio power that 
were more pronounced at smaller spatial scales. At these 
smaller scales they found that the correlation amplitude in- 
creased with radio power up to ~ 10 25 W/Hz, above which 
the correlation amplitude flattened and then began to fall 
again above ~ 10 25 ' 5 W/Hz. 

We divide our sample of radio LRGs into three bins 
of roughly equal size by their flux density with cuts at 
Si.4GHz =5 and 20mJy/Hz. Fig. [3] shows the LxR cross 
correlations for each of the flux density bins for the three 
LRG samples. 

We do not find the same behaviour in each LRG sam- 
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Figure 2. The LxL and LxR correlation functions for the SDSS, 2SLAQ and AAS1 samples (left to right). Also shown is the LxM cross 
correlation of the sample matched to the radio samples in terms of optical luminosity and colour. The top panels show the correlation 
functions along with best fit models to the LxM (solid) and LxR (dashed) data. The models are double power laws in £(r) mapped to 
w(8) with Limber's equation. The bottom panels show the ratio of the LxR and LxM correlations to the best fit LxL model. Errorbars 
are omitted from the top plots for clarity. 
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Figure 3. LxR cross correlation functions for the SDSS, 2SLAQ and AAC samples (left to right) with the radio sample split into 
three bins by radio flux density. The bottom panels show the ratio of the correlation functions to the best fit LxL model. Errorbars are 
omitted from the top plots for clarity. In the SDSS sample there is evidence for a positive correlation between radio flux density and the 
correlation strength. This is not so clear in the other two samples. 



pie. In the SDSS sample there is a clear positive correlation 
between radio luminosity and the angular correlation am- 
plitude. As with the difference between the LxL and LxR 
correlations, the luminosity dependence is most pronounced 
at small angular scales (<2 arcmin). In the 2SLAQ sam- 
ple the lowest flux density bin is significantly (> 3a) lower 
than the other two subsets that essentially lie on top of each 
other. However, in the highest redshift AAJ1 sample we find 
no significant difference between the three flux bins. 



Our resu l ts ar e broadly consistent with those of 
iDonoso et al.l |2009). At the median redshift of our sam- 
ples the flux density limits we have chosen 5/20 mjy corre- 
spond to log(Li.4GHz/(W/Hz)) = 24.3/24.9, 24.7/25.3 and 
25.0/25.5 for the SDSS, 2SLAQ and AAfi samples respec- 
tively. Hence our results are consistent with the small-scale 
radio correlation function increasing for radio luminosities 
< 10 25 W/Hz and then flattening off. 
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Figure 4. The SDSS* x SDSS", 2SLAQ* X2SLAQ" and 
AAf2*xAAQ" cross correlation functions. The solid lines show 
the best-fit double power laws in f(r) mapped to w(8) with Lim- 
ber's equation. 



6 CLUSTERING STRENGTH OF RADIO 
SOURCES AT Z < 0.68 

Given our three samples span z ~ 0.68 to 0.35 (~2.3 Gyr of 
cosmic time) we are interested in any signs of evolution in 
the clustering strength of radio LRGs. To make a fair com- 
parison we make our samples as equivale nt as possible in 
terms of their optical and radio selection. ISawangwit et al.l 
(2009) showed that by applying nragnitude limits of idev = 
19.32 and 20.25 to the 2SLAQ and AA£2 samples respec- 
tively, they could be made roughly equivalent to the SDSS 
sample in terms of their (optical) luminosity distribution 
and space density (given that the LRG luminosi ty function 
does not evolve strongly over this interval; e.g. IWake et al.l 
120081 ). They call these optical LRG samples with additional 
magnitude cuts SDSS*, 2SLAQ* and AAf2* . In addition, we 
increase the radio flux limit of the SDSS* and 2SLAQ* sam- 
ples to 13.8 and 4.9 mjy respectively, in order to be equiva- 
lent to a radio power of ~ 10 24 ' 72 W/Hz at the median red- 
shift of each sample. This results in 3576, 4755, 5089 objects 
in our radio matched sanrples we will call SDSS", 2SLAQ" 
and AAfi" . 

In Fig. [4] we show the cross correlations between the 
radio nratched and optically matched samples along with a 
best-fit double power law in £(r) mapped to w(0) with Lim- 
ber's equation. Again the break in the power law is assumed 
to be at l/i" 1 Mpc when performing the fit, in addition to 
which we fix the slope of the large-scale power law to be 

7 = 1.8 to facilitate comparisons of the clustering ampli- 
tude. 

Fig. \5\ shows the large-scale (> l/i -1 Mpc) clustering 
amplitude as a function of redshift (squares). The results are 
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Figure 5. The large-scale cross-correlation amplitude for the 
SDSS*xSDSS", 2SLAQ*x2SLAQ" and AAfTxAAQ" corre- 
lations (squares) and the SDSS* x SDSS 26 *, 2SLAQ* X2SLAQ 26 * 
and AAQ*xAAf2 26 * correlations (circles; these have been offset 
horizontally for clarity). Neither show any evidence for a trend 
with redshift. The lines show models for the cross-correlation am- 
plitude based on differing assumptions about the radio-source au- 
tocorrelation. The solid and dashed lines are constant dark halo 
mass and long lived models respectively, fit to the lower luminos- 
ity (square) points. The dot-dashed line assumes that the radio 
autocorrelation is the same as that measured for quasars. 



consistent with a constant clustering amplitude with little 
evidence for evolution. 

To make a comparison with previous work on the clus- 
tering of radio sources we calculate the radio LRG two-point 
autocorrelation function in each sample and show the results 
in Fig. [6] The lines in the figure show the best fit power law 
in £(r) (7 fixed to 1.8) mapped to w(9) with Limber's equa- 
tion. Table [2] summarises the results of our power law fitting 
in t his section. 

iBrand et al.l (|2005l ) measured the clustering strength 
of low- luminosity radio galaxies at 2 ~ 0.3 and found 
r = 6.1 ± 1.1 ft -1 Mpc assuming £(r) = (r/ro) -18 . This 
is significantly (~ 3.1 a) less than the value of rp = 11.0 ± 
Peacock fc Nicholson! (|l99ll ) at 



z < 



1.2 ft" 1 Mpc found by 

0.1. IBrand et al.l ((2005) explained some of the discrepancy 



by assuming the clustering amplitude would increase accord - 
ingly with linear theory from z = 0.3 to O. lWake et al.l (120081 ) 
derived an autocorrelation clustering strength (extrapolated 
from the cross correlation) of ro = 12. 3± 1.2 hT 1 Mpc for ra- 
dio sources in the 2SLAQ sample. 

In Fig [7] we show the values for ro calculated 
in this paper (circles) alo ng with t he values found by 
iPeacock fc Nicholson! (ll99lD (tr iangle). IBrand et all (120051 ) 
(square) and IWake et al.l ( 20081 ) (star). A best fit to these 
points gives a gradient co nsistent with zero (< . 5cr). 

The flux limit s of the IPeacock fc Nicholson! (|l99ll ) and 
(jWake et al.ll2008l ) samples are such that they sample < 2 
times fainter in terms of radio luminosity at their median 
redshifts compared to our samples. Hence, these surveys are 
sampling relatively bright objects where we do not find clus- 
tering strength cor relating with radio power (Fig. [3}. The 
IBrand et al.l (|2005l ) sample on the other hand goes ~ 6.5 
times fainter. This difference may explain why their result 
is lower than the other values. 
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Table 2. Parameters of power law fits to the angular correlation function for samples matched in terms of their optical and radio 
properties. Samples included are denoted * (optically selected LRGs), r* (radio LRGs with Li.4GHz > 10 W/Hz) and 26* (radio 

LRGs with Li.4GHz > 10 26 W/Hz). 



Correlation 


SDSS 
r (h~ 1 Mpc) 


7 


2SLAQ 
r ( )(/i _1 Mpc) 


7 


AAH 
r (h~ 1 Mpc) 


7 


* x r* 


4.55 
10.67±0.28 


2.81±0.09 
1.8 


5.07 2.72±0.08 
11.63±0.45 1.8 


5.49 
10.96±0.44 


2.53±0.06 
1.8 


V * XT'* 


10.4±1.6 


1.8 


11.9±1.3 


1.8 


11.6±1.3 


1.8 


* X 26* 


10.67±1.13 


1.8 


11.63±0.99 


1.8 


10.96±0.79 


1.8 
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Figure 6. Autocorrelation functions for the three samples after 
they have been made approximately equivalent in terms of their 
optical and radio selection. The solid line gives the best fit relation 
in terms of £(r) = (ro/r) 7 propagated through Limber's formula. 



6.1 Evolution of the clustering amplitude of radio 
galaxies 

It is clear from Fig. [7] and the squares in Fig. [5] that there 
is little evidence for evolution in the clustering amplitude 
of either the auto or cross correlation with redshift. Fitting 
ro oc (1 + z) a to the measurements in Fig. [7] (excluding 
the iBrand et ail 120051 point) gives a = 0.20 ± 0.30. The 
same fit for the cross correlation (squares in Fig. [S| gives 
a — 0.23 ± 0.21. In addition we fit a series of physically 
motivated models to our data: 

(i) In linear theory structures evolve as £ m ex D 2 (z) where 
D(z) is the growth factor and £ m is the matter correlation 
function. Structures that grow according to this have a con- 
stant bias (6) where b 2 = £/£ m - In Fig. [7] we fit a constant 
bias model to our data and plot the results as the dot-dashed 
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Figure 7. The radio galaxy X radio galaxy autocorrelation ampli- 
tude ro as a function of redshift. The circ les correspond to our 
three samples in addition to points from iPeacock &: Nicholson 



three samples in addition to points from Fcacock tk Nicholsor 
1 199lh f triangle), brand et ait J2005h (square) and IWake et~ 



12008) (star; this point is offset horizontally for clarity). The 
measurements are consistent with no evolution in the clustcr- 
ing am plitude with redshift up to z ~ 0.68. The IBrand et al.l 
(2005) point is significantly offset from the rest, potentially due 
to their sample being dominated by fainter o bjects. The lines give 
three model fits to the data excluding the IBrand et al,l (2005) 
point. The models are constant dark halo mass (solid), long-lived 
(dashed) and linear growth (dot-dashed). 



line. This model is only marginally consistent with our data 
( X 2 = ll.ll;P(xV = 4)~2.5%). 

(ii) We further include the long-lived model of lFrvl (119961 ) 
that assumes no change in the comoving number density 
of galaxies and that clustering grows solely due to gravi- 
tational effects. The long lived model has been shown to 
describe evolution in t he LRGxLRG correlation function 
(jSawangwit et alj|2009l ). Fits to the radio cross and auto 
correlation amplitudes are shown as dashed lines in Figs. [5] 
and [7] Note that the long-lived model does not require that 
radio sources themselves exist for cosmological timescales. 
Rather that the dark matter halos that host them evolve 
passively under the influence of gravity. This model involves 
little evolution and hence gives a better fit to our data. 

(iii) We also fit a model in which dark halo mass is 
held constant with re dshift (solid lines Figs. [5] and [7] see 
ISawangwit et al.ll2009l for our formalism for relating bias to 
dark halo mass). This type of model has been shown to fit 
the evolution of the quasar autocorrelation function over 
< z < 2.5. 
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To relate our models for the radio auto correlation to 
the cross correlation amplitudes in Fig. [5] we assume the 
LRG clustering f ollows the best-fit long-lived model from 
ISawangwit et al.l ((2009J). We then assume that the cross- 
correlation of two samples (£12) is related to the tw o au- 
tocorrelations (£n and £22) by (e.g. IWake et al.ll2008l ) 



£12 = £1162 



(2) 



The long-lived and constant mass models fit our data 
acceptably {P(x 2 , v) > 10 %) with the exception of the long- 
lived model fit to the cross correlation results. In this case 
we find x* = 9.08 with P(x 2 ,^ = 2) ~ 1.1% indicating 
that the long-lived model may not accurately describe radio 
source clustering. The dark halo masses derived for the radio 
sources are 9.4 and 5.6 x 10 13 /t -1 M© for the cross and auto 
correlation respectively. 

The offset in the derived dark halo masses may be due 
to equation [2] no t hold ing in our data. In this paper and 
ISawangwit et al.l (|2009r ) large-scale correlation amplitudes 
have been calculated for the * x * LRG auto correlation, 
* x *r cross correlation and the *r x *r autocorrelation. Cal- 
culating the cross correlation amplitude from the autocorre- 
lations and equation[2]gives ro — 9.75±0.75, 10.28±0.56 and 
10.26iO.58/iT 1 Mpc for the SDSS, 2SLAQ and AAfi samples 
respectively. These values are consistently lower than the 
amplitudes measured from the cross correlation, although 
only by ~ 1 to 2a. Taking the three samples together the 
offset is ~ 3(J. The offset may be due to two causes. First, 
the fits are not completely equivalent; we fit double power 
laws over a larger range of 9 for the * x * and * x *r correla- 
tions. Second, equation [2] assumes a linear bias model that 
may not hold in this case. 
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Figure 8. The SDSS* x SDSS 26 *, 2SLAQ* X2SLAQ 26 * and 
AAQ*xAAf2 26 * cross correlation functions. The solid line shows 
the best fit solutions to £(r) = (ro/r) 18 propagated through Lim- 
ber's formula. 



I Johnston et al.l 120081 ). Hence, while an LRG selection finds 
most hosts of radio galaxies in the range 10 23 < L < 
10 26 W/Hz, above this there is a significant population of 
blue objects that may be missed by the optical selection. 



7 CLUSTERING OF LUMINOUS RADIO 
SOURCES AT Z < 0.68 

The SDSS", 2SLAQ" and AA£3 r * samples are dominated 
by objects below the ~ 10 26 W/Hz divide between brighter 
FRI and FRII sources. 10 26 W/Hz is also roughly the divide 
between radio AGN host galaxies that show no emission 
lines and those that do. Furthermore, above ~ 10 26 W/Hz 
radio sources show strong number density evolution, more 
similar to quasars. Th e clustering strength of quasars evolve s 
weakly with redshift (jCroom et al.ll2005l : IRoss et al.ll2009h . 
In this section we use our sample to test for evolution in the 
clustering amplitude of bright radio LRGs. 

A radio luminosity of 10 26 W/Hz corresponds to flux 
limits of ~ 266, 95 and 57mJy at z ~ 0.35, 0.55 and 
0.68 respectively. Cutting back our samples to these lim- 
its leave 49, 224 and 283 objects in what we will call the 
SDSS 26 *, 2SLAQ 26 * and AAfi 26 * samples. These are not 
large enough to calculate an autocorrelation amplitude so 
we calculate the two-point cross correlation with the corre- 
sponding (starred) LRG catalogues. Fig. [8] shows the cross 
correlations along with best fit solutions to f (r) = (ro/r) 1,8 
propagated through Limber's formula. In Fig. [S] the circles 
show ro for these bright objects as a function of redshift. 

Note that while lower-luminosity radio AGN are typi- 
cally found in LRGs with red stellar spectra, the brightest 
sources can be found in both broad and narrow emission- 
line objects that have younger stellar populations (e.g. 



7.1 Evolution of the clustering amplitude of 
luminous radio sources 

As with the less luminous sample, we find very little evidence 
for evolution in the clustering amplitude of these brightest 
radio sources. Fitting the cross correlation points (circles 
Fig. 0) for r oc (1 + z) a gives a = -0.38 ± 0.52. 

We also compare our results for the clustering of the 
brightest radio sources with that of quasars. To do this we 
assume all quasars cluster such that the ir derived dark halo 
mass is 3 x 10 1 2 M at all redshifts (|Croom et all 120051 ; 
IRoss et all 120091 ) and work backwards to derive the clus- 
tering strength as a function of redshift. We then convert 
this to a cross-correlation clustering strength making the 
assumptions discussed above. We plot the derived value for 
the cross-correlation amplitude in Fig. [5] as the dot-dashed 
line. 

While our derivation is approximate it suggests that 
the clustering amplitude of bright radio sources may 
evolve slowly with redshift and in an almost luminosity- 
independent manner, similar to the results found for 
quasars. There is also strong evidence that the objects in our 
sample cluster considerably more strongly than quasars, as 
previously noted for fainter radio sources. The implication 
is that these bright radio sources must inhabit considerably 
denser regions of space and correspondingly more massive 
dark matter halos compared to quasars. There is less differ- 
ence between the large-scale clustering environment of lu- 
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minous and faint radio sources. However, the small scale 
r < lh _1 Mpc clustering environment of the brightest radio 
sources appears significantly denser than the environment 
of the faintest sources. This may indicate the effect of the 
radio source on its immediate surroundings, given that the 
larger-scale environment is luminosity independent. 



8 CONCLUSIONS 

We have cross-matched three large photometric samples of 
LRGs to the FIRST and NVSS surveys to define three sam- 
ples of radio LRGs. We have then measured the 2-point cor- 
relation function to investigate the clustering properties of 
the samples. Similar to previous authors we find evidence 
that: 1) radio LRGs are more strongly clustered than non- 
radio galaxies matched in terms of luminosity and colour; 2) 
the clustering strength of radio LRGs at r < Ih _1 Mpc in- 
creases with radio luminosity up to L ~ 10 25 W/Hz and then 
remains roughly constant but at r > lh _1 Mpc the cluster- 
ing is independent of radio luminosity. We further find that 
inconsistencies between the radio-LRG cross-correlation and 
auto-correlation amplitudes may suggest that a simple lin- 
ear bias model may be insufficient to describe the relation 
between radio-LRG and/or LRG clustering and mass clus- 
tering. 

However, the primary goal of this work was to investi- 
gate any evolution in the clustering strength of radio galaxies 
with redshift. We find no evidence for evolution in the large 
scale (r > lh _1 Mpc) clustering amplitude of radio-LRGs. 
We show that our radio x radio LRG autocorrelations are 
consistent with previous authors indicating no evolution in 
comoving coordinates in the redshift range < z < 0.68 
(~ 6Gyr). Furthermore, we make use of cross-correlations 
to increase the precision of our results and still find no ev- 
idence for evolution. We then restrict our samples to ob- 
jects with L > 10 26 W/Hz and again find no evidence for 
evolution in the correlation amplitude. Our results are con- 
sistent with a single dark halo mass of 9 x 10 13 /i _1 Mq for 
all radio LRGs in the redshift range < z < 0.68. This 
could be compared to quasars that appear to inhabit halos 
of uidh = 3x 10 12 /i _1 Mq again reasonably independent of 
redshift and luminosity. The significantly more massive halos 
for the high-luminosity radio galaxies, at least, may provide 
a problem for unified AGN models. Finally, our radio-LRG 
cross-correlations are inconsistent with a model in which 
LRG clustering follows a long-lived model and radio sources 
are randomly sampling the LRGs. 



Brown M. J. I., et al., 2007, ApJ, 654, 858 
Cannon R., et al., 2006, MNRAS, 372, 425 
Clewley L., Jarvis M. J., 2004, MNRAS, 352, 909 
Condon J. J., Cotton W. D., Greisen E. W., Yin Q. F., 

Perley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 

1693 
Croom S. M., et al., 2005, MNRAS, 356, 415 
Croom S. M., et al., 2009, MNRAS, 399, 1755 
Donoso E., Best P. N., Kauffmann G., 2009, MNRAS, 392, 

617 
Dunlop J. S., Peacock J. A., 1990, MNRAS, 247, 19 
Eisenstein D. J., et al., 2001, AJ, 122, 2267 
Fanaroff B. L., Riley J. M., 1974, MNRAS, 167, 31P 
Fry J. N., 1996, ApJL, 461, L65+ 
Hamilton A. J. S., 1993, ApJ, 417, 19 
Hickox R. C, et al., 2009, ApJ, 696, 891 
Hine R. G., Longair M. S., 1979, MNRAS, 188, 111 
Johnston H. M., Sadler E. M., Cannon R., Croom S. M., 

Ross N. P., Schneider D. P., 2008, MNRAS, 384, 692 
Kimball A. E., Ivezic Z., 2008, AJ, 136, 684 
Landy S. D., Szalay A. S., 1993, ApJ, 412, 64 
Longair M. S., 1966, MNRAS, 133, 421 
Magliocchetti M., et al., 2004, MNRAS, 350, 1485 
Mandelbaum R., Li C, Kauffmann G., White S. D. M., 

2009, MNRAS, 393, 377 
Mauch T., Sadler E. M., 2007, MNRAS, 375, 931 
McCarthy P. J., 1993, ARAA, 31, 639 
Peacock J. A., Nicholson D., 1991, MNRAS, 253, 307 
Phillipps S., Fong R., Fall R. S. E. S. M., MacGillivray 

H. T., 1978, MNRAS, 182, 673 
Prestage R. M., Peacock J. A., 1988, MNRAS, 230, 131 
Ross N. P., Shanks T., Cannon R. D., Wake D. A., Sharp 

R. G., Croom S. M., Peacock J. A., 2008, MNRAS, 387, 

1323 
Ross N. P., et al., 2009, ApJ, 697, 1634 
Sadler E. M., et al., 2007, MNRAS, 381, 211 
Sawangwit U., et al., 2009, ArXiv e-prints 
Schmidt M., 1968, ApJ, 151, 393 
Wake D. A., Croom S. M., Sadler E. M., Johnston H. M., 

2008, MNRAS, 391, 1674 
Wake D. A., et al., 2006, MNRAS, 372, 537 
Wake D. A., et al., 2008, MNRAS, 387, 1045 
Yates M. G., Miller L., Peacock J. A., 1989, MNRAS, 240, 

129 
York D. G., et al., 2000, AJ, 120, 1579 



REFERENCES 

Abazajian K. N., et al., 2009, ApJS, 182, 543 

Adelman-McCarthy J. K., et al., 2007, ApJS, 172, 634 

Becker R. H., White R. L., Gregg M. D., Brotherton M. S., 
Laurent-Muehleisen S. A., Arav N., 2000, ApJ, 538, 72 

Bell E. F., et al., 2004, ApJ, 608, 752 

Best P. N., Kauffmann G., Heckman T. M., Ivezic Z., 2005, 
MNRAS, 362, 9 

Blake C, Wall J., 2002, MNRAS, 329, L37 

Brand K., Rawlings S., Hill G. J., Tufts J. R., 2005, MN- 
RAS, 357, 1231 



